SCIENTIFIC 

REPORTS 



s 




SUBJECT AREAS: 

CHEMICAL PHYSICS 

COMPUTATIONAL CHEMISTRY 

CARBON NANOTUBES AND 
FULLERENES 

NANOPARTICLES 

Received 
8 January 2013 

Accepted 
12 February 2013 

Published 
26 February 201 3 



Correspondence and 
requests for materials 
should be addressed to 
P.M. (huaiping@sinap. 
ac.cn); R.Q.Z. 
(a p rqz@c i ty u . ed u . h k) 
orZ.G.W. (wangzg@ 
ilu. edu.cn) 



Defect Induced Electronic Structure of 
Uranofullerene 

Xing Da\\ Cheng Cheng^, Wei Zhang^, Minsi Xin^ "^, Ping Huai^, Ruiqin Zhang^"^ & Zhigang Wang^'^ 

^ Institute of Atomic and Molecular Physicsjilin University, Changchun, 1 3001 2, P. R. China, ^Shanghai Institute of Applied Physics, 
Chinese Academy of Sciences, Shanghai 201 800, P. R. China, ^Department of Physics and Materials Science and Centre for 
Functional Photonics (CFP), City University of Hong Kong, Hong Kong SAR, P. R. China, ^Beijing Computational Science Research 
Center, Beijing, P. R. China, ^State Key Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry 
Jilin University, Changchun, 130023, P. R. China. 

The interaction between the inner atoms/cluster and the outer fuUerene cage is the source of various novel 
properties of endohedral metallofuUerenes. Herein, we introduce an adatom-type spin polarization defect 
on the surface of a typical endohedral stable U2@C6o to predict the associated structure and electronic 
properties of U2@C6i based on the density functional theory method. We found that defect induces obvious 
changes in the electronic structure of this metallofuUerene. More interestingly, the ground state of U2@C6i is 
nonet spin in contrast to the septet of U2@C6o. Electronic structure analysis shows that the inner U atoms 
and the C ad-atom on the surface of the cage contribute together to this spin state, which is brought about by 
a ferromagnetic coupling between the spin of the unpaired electrons of the U atoms and the C ad-atom. This 
discovery may provide a possible approach to adapt the electronic structure properties of endohedral 
metallofuUerenes. 

Endohedral metallofuUerenes (EMFs) are novel fullerene derivatives that can encapsulate metal atom or even 
clusters in their inner space^'^. Previous studies show that the inner cluster and/or the outer cage of the 
unique structure preferentially adopt the stable closed-shell electronic configuration^'^, such as La2@C72^, 
La2@C8o'', Sc2@C66', Sc2C2@C68', Sc3N@Cn (n = 68, 78, 80)^'-^^ Ti2C2@C78'^ There is a significant charge 
transfer between the metal atoms and the fullerene In EMFs, the internal metal atoms always donate electrons 
to the fullerene cage and carry considerable positive charges. This can result in an enhanced chemical reactivity as 
compared to the hollow fullerene^'^. EMFs have strong ability of donate and accept electrons, implying that they 
can be easily oxidized and reduced^^. The stability of the whole system can be altered by this oxidation or 
reduction. The physical properties of EMFs strongly depend on the properties of the inner metal atom. Such a 
series of fantastic M@Cn structures have drawn extensive research interests in recent years. It is well known that 
defects can be easily formed in a large variety of carbon nanostructures, such as fullerenes, nanotubes and 
graphene. This defect can also bring about a significant change in the electronic structure as well as many other 
properties^^"^°. But the natures of defective EMFs system are relatively unexplored till date. This work is aimed at 
gaining an understanding in this area through the study of an EMF structure with a typical spin polarized 
defective fullerene, which will be potentially beneficial in the design or synthesis of novel materials in several 
areas including catalysis and medicine^ 

EMFs contain U, Th, Pu and other actinide elements are very important complexes. Actinide metallofuUerene 
materials not only have the thermal stability of the fullerene-like, but also have the physicochemical properties of 
the actinide elements, and therefore it is foreseeable that they will have important potential application in 
advanced fuel technology. So, studies on EMFs containing U atom have been widely reported. Guo et al. had 
prepared uranium fullerenes U@C28 and U@C6o experimentally for the first time. C28 is usually unstable. The U 
atom wrapped into C28 stabilizes the fullerene cage^^. Subsequently, Diener et al. synthesized U@Cy4 and U@ 
C82^^- Zhao and his co-workers successfully synthesized and isolated a new type of metallofuUerene containing 
metal atoms of transuranic elements for the first time, proving that actinide metal atoms can stabilize the 
otherwise unstable C82 fullerene^^ and explained the structural stability of U@C82 through theoretical calculation 
of density functional theory (DFT)^^. The uranium dimer compound is a rare compound whose presence is 
evidenced only in experimental spectra^^. Therefore, its theoretical study is scarce. Gagliardi and Roos have 
theoretically studied the U2 molecule using high-level ab initio method and found that this molecule contains 
three electron-pair bonds and four one-electron bonds, corresponding to a unique quintuple bond, and the U2 
molecule has a septet ground state^^. Theorists have also predicted the possibility of fullerene encapsulating two U 
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atoms based on quantum first principles method. Usually, the exo- 
hedral dimetal fullerenes are more stable than the endohedral struc- 
tures, e.g. Cr2, M02 and W2^^. However, the U2 case is different. Its 
endohedral complex is significantly more stable than the exohedral 
and U2 lies at the center of the C^q fullerene^^. The U-U metal bond 
tends to dissociate and the U atoms interact with the inner wall of the 
cage when the neutral U2 molecule is encapsulated into carbon cage; 
but due to the special size of the Ceo fullerene, the U-U bond is 
weakened but not completely broken in the U2@C6o EMF. The 
results from a DFT study using the GGA-PBE exchange correlation 
functional shows that the dimetalloendofullerene, U2@C6o, has an 
unprecedented U-U multiple bond consisting solely of sixfold ferro- 
magnetically coupled one-electron-two-center bonds, which are 
dominated by the uranium 5f atomic orbitals, and it is septet in 
ground state^°. This work opened up a new chemical field of the 
interaction between the actinide metal-metal bond and the fullerene, 
and provided the basis for many other studies. The size of the cavity 
can affect the properties of EMFs. If we enlarge the size to, say, C70 
and C84, the metal-metal bond will be completely dissociated but the 
whole systems still have a septet ground state^^. 

Carbon nanostructures occurring in nature mostly contain 
defects, which cause changes in the electronic structure of carbon 
based nanomaterials and have been widely reported. However, full- 
erene is a necessary ingredient in EMFs. If the cage surface of the 
EMF which encapsulates U2 molecule within it contains a defect, the 
nature of such a structure is expected to be much different. So, we 
have considered an adatom defect by introducing a carbon ad- atom 
over the [6, 6] bond on the typical U2@C6o system and analyzed the 
features of this defective EMF U2@C6i. 

Results 

In order to examine the reliability of the exchange correlation (XC) 
functional and basis sets, we calculated the isolated U2 molecule, the 
Cgo fullerene and the EMF, U2@C6o- All the results qualitatively agree 
with the previous reports. Furthermore, the results showed that 
defective C^i has a triplet ground state. The adatom contributes 
significantly to the spin. The details for the calculations have been 
presented in the supplementary information Part 1.3. 

In Cei fullerene, we have considered six isomers with the two 
uranium atoms at different positions. The full geometry optimization 
results showed that the total energy of the whole system is progres- 
sively reduced as a uranium atom gradually moves closer to the 
defect. Eventually, we found the structure c in all isomers to be the 
most stable. Structure c can be characterized as one U atom located 
below the defect and another U atom facing a six-ring as shown in 
Figure 1. The details for searching the lowest energy structure are 
provided in the supplementary information Part 2. 
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Figure 1 | The equilibrium structure (Structure c) with the lowest energy 
of all the U2@C6i isomers. 



To investigate the possible mechanism responsible for the struc- 
tural stability of c, we have proceeded systematically. First of all, a 
qualitative explanation is provided as follows. The single U atom is 
reported to tend to adsorb on a hexagonal face from the inside when 
it is encapsulated in a carbon cage^^. The neutral U2 molecule is found 
to show a tendency to dissociate and also adsorb on the hexagonal 
faces in C^o fullerene^^'^^. Interestingly, the adatom-type defect can 
enlarge the local size of the C^q. We have calculated and compared 
the size of the C^i and C^q. The diameter of C^i is about 0.31 A longer 
than the original C^q (Ih) along the direction of the defect. Moreover, 
the defect dissociates the original [6, 6] bond of the C^q. 
Consequently, the distance between these two C atoms formerly 
forming the [6, 6] bond in the original Ceo increases by 0.79 A. 
This causes the inner U atom to sense the outer unsaturated ad- atom. 
Based on the discussion above, there is no doubt that both the size 
effect and the impact of the adatom directionally guide the motion of 
one of the two U atoms to the defect area. As a result, the whole 
system is energetically favored to form the unique structure c as 
depicted in Figure 1. The distance of the two U atoms in C^i fullerene 
(2.80 A of BP86, 2.88 A of PBEO, see Table 1) is larger than the U2@ 
Ceo case which has been calculated at the same level (2.70 A of BP86, 
see supplementary information Part 1.3), including all the U2@C6o 
results reported earlier (2.72 A of PBE'°, 2.73 A of BP86''). This also 
confirms the conclusions of Gagliardi etc.^^, that is, if the space 
around U2 in U2@Cn becomes larger, the two U atoms will tend to 
dissociate and interact with the inner wall of the cage^^. It should be 
noted that, one of the adsorption sites of the two U atoms in Cei 
fullerene is different from that of U2@C6o, the C atoms in the defect- 
ive area form a 10-membered ring-like structure and a U atom shifts 
below the center of this 10-membered ring. This position is different 
from the case of U2@C6o in which the two U atoms face the centers of 
two hexagonal rings. In this way, we have arrived at an interesting 
'dissociation- adsorption' equilibrium state structure of neutral U2 
molecule encapsulated in defective C^i fullerene. In Part 1.3 of the 
supplementary information, we have explained all the details of size 
comparison between C^q and C^i fullerenes. 

Table 1 shows the full geometry optimization results of U2@C6i 
(structure c) calculated at BP86/ANO_ECP60MWB~3-21G level 
and also with the PBEO data as a test. From the energy relations of 
different electronic state we can see both the two XC functional 
show that the nonet state is the most stable state and there is no 
imaginary frequency in the frequency calculation. Hence, structure 
c is identified as the ground state structure and we express this 
electronic state as U2@C6i (c, ^A'). 

Since infrared (IR) and Raman spectroscopy have been widely 
used to characterize metallofullerenesl, we showed the IR and 
Raman spectrums of U2@C6i (c, ^A') predicted at BP86/ 
ANO_ECP60MWB~3-21G level in this paper. These results are 
expected to guide experimental research in this area. We first exam- 
ined the infrared and Raman spectrums of the C^q and U2@C6o and 
the results are consistent with the previous studies (see supplement- 
ary information Part 1.2). Also, we performed the same calculation 
for defective C^i fullerene. The IR-active vibrational mode in Cei is 
more complex than that in C^q. Unique peaks due to the ad- atom 
appear in the IR spectrum at the frequencies of 1257.5 and 



Table 1 


Full geometry optimization 


results of U2@C6i 


(structure c) 


System 


Method 


El. state 


Total energy (H) 


Ru-u (A) 


U2@C6i 


(c) BP86 


^A' 


-3265.9161 1 


2.77 








-3265.91672 


2.80 








-3265.89768 


2.79 




PBEO 


'A" 


-3262.79316 


2.84 






^A' 


-3262.79947 


2.88 






11 A. 


-3262.78450 


2.71 
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Figure 2 | Simulated IR and Raman spectrums of U2@C6i (c, ^A') at BP86/ANO_ECP60MWB~3-21G level. 



1353.6 cm~\ These vibrational modes correspond to the asymmetric 
and symmetric stretching modes of the two C-C bonds of the ad- 
atom respectively (see supplementary information Part 1.2). 

Figure 2 shows the IR and Raman spectra of U2@C6i (c, ^A') 
respectively which are predicted at the BP86/ANO_ECP60MWB~ 
3-21G level. We examined all the vibrational mode related to the U 
atoms and found that the symmetric U-U stretching mode is IR- 
inactive but Raman-active. So, we can find this signal in the 
Raman spectrum at the frequency of 127.4 cm~\ This mode can 
suggest that there is still interaction between the two U atoms in 
Cei fullerene, the special space size of C^i may not cause them to 
dissociate completely, although the distance between the two U 
atoms was enlarged in C^i. The asymmetric U-U stretching mode 
is IR- active and its weak intensity is predicted at a frequency of 
142.0 cm"\ This mode explains the interaction between the U atoms 
and the inner wall of the cage. The unique IR-active peaks arising due 
to the adatom lie at the frequencies of 1327.5 and 1341.6 cm"\ The 
vibrational modes still correspond to the asymmetric and symmetric 
stretching modes of the two C-C bond of the adatom respectively. 

Furthermore, we compared the vibrational frequencies of U atoms 
between U2@C6o and U2@C6i systems. From Table 2 we can clearly 
see that the relative frequency positions in the two EMF systems are 
reversed. The frequency of symmetric U-U stretching mode in U2@ 
Cei becomes significantly smaller than that in the U2@C6o. This 



Table 2 | The frequencies of symmetric and asymmetric U-U 
stretching modes in U2@C6o and U2@C6i systems. All of these 
results were calculated at BP86/ANO_ECP60MWB~3-21 G level 



System 



Symmetric U-U stretching Asymmetric U-U stretching 



U2@C60 
U2@C61 (C, 



A') 



162.9 
127.4 



157.2 
142.0 



interesting phenomenon illustrate that the bond between the two 
U atoms in C^i is very weak. 

The key factor that governs the stability of EMFs is the magnitude 
of electronic charge transfer between the inner atom/cluster and the 
cage. Therefore, it is necessary to analyze its charge distribution. We 
have named the adatom as C^uh^ the U atom close to the defect as 
Uist, and the other U atom as U2nd- We have particularly focused on 
the charge distribution on these three atoms and the relevant data is 
listed in Table 3. Usually, metal atoms in fullerene, including the U 
atom, donate electrons to the cage and attain a considerable positive 
charge. Owing to the different equilibrium positions of the two U 
atoms in the EMF U2@C6i, Uist and U2nd also carry different amount 
of positive charges. In Table 3, we observe that the Uist is more 
positively charged than the U2nd- It means that Uist which lies below 
the defect area transfers more electronic charge to the fullerene cage 
than the U2nd which faces a six-member ring in the C^i fullerene. 
This result can be explained as follows. The Uist adsorbs on a 10- 
membered ring-like structure in contrast to the U2nd which faces a 
hexagonal structure. Uist interacts with more C atoms than the U2nd- 
As a result, it donates more electronic charge. The C atoms of the 
cage gather electronic charge from the inner U atoms and thus get 
negatively charged except the adatom. Both the two XC-functionals 
show that the C ad-atom transfers a small amount of 



Table 3 Mulliken charge distribution 


of U2@ 


C61 (c,'A') 


System Method 


Atom 


Charge (e) 


U2@C6i (c, ^A') BP86 


C61st 


0.122 




Uu, 


2.383 




U2nd 


2.316 


PBEO 


C61st 


0.181 




U,s, 


2.491 




U2„d 


2.410 
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0.500 
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0.357 



0.286 



0.214 



0.143 



0.071 



0.000 



-5.51 



-3.91 



-2.31 -0.71 0.89 2.49 
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4.09 5.69 



Figure 3 | The electron density color-filled map of U2@C6i (c, ^A'). 



electronic charge to the cage and acquires positive charge like the 
inner U atoms. The ten C atoms of the 10-membered ring mentioned 
above are comparatively more negatively charged (about —0.297) 
than other C atoms on the cage, as they have gained electronic charge 
from both the inner U atom and the outer C adatom. 

Electron density map shows the charge distribution in the system 
more clearly. In Figure 3, the region between the U atoms and the 
inner wall of the cage shows apparent electron density, suggesting 
obvious bonding interaction in these areas. The electron density in 
the space between the Uist and Ceist is very low. We have also given 



the electron density contour line map of U2@C6i (c, ^A') as Figure S6 
in the supplementary information in Part 1.3. 

Molecular orbital (MO) analysis also helps us to understand the 
electronic structure. The multiplicity of the ground state of U2@C6i is 
M = 9, which corresponds to eight unpaired electrons. So, we are 
concerned with the characteristics of LUMOa,HOMOot,HOMOot-l 
to HOMOa-7, which are singly- occupied frontier molecular orbitals 
of the U2@C6i(c, ^A'). It can be seen from Figure 4 that the two U 
atoms provide significant contribution to all these orbitals. Due to 
impact of adatom-type defect on the fullerene, the shape of the 
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Table 4 | The percentages of atomic orbitals active in intermoleculcr interaction 

Uut U2„d 



U9@Qi (c, ^A') 


5f 


6d 


5f 


6d 


2s 


2p 




DU.DV /o 


n mo/ 
U.Uo /o 


O C Q CO/ 

zO.oD /o 


n 1 AO/ 

U . \0/o 


n mo/ 
U.Uo /o 


n ^"70/ 
U.D/ /o 


HOMOoc 


24.97% 


0.63% 


46.62% 


0.24% 


0.09% 


1 .09% 


HOMOoc-1 


40.01% 


1.17% 


22.54% 


0.18% 


1.15% 


6.30% 


HOMOoc-2 


29.94% 


0.31% 


40.72% 


0.23% 


0 


0 


HOMOoc-3 


19.52% 


0.52% 


26.16% 


0.02% 


0.65% 


2.77% 


HOMOoc-4 


12.84% 


0.01% 


21.55% 


0.06% 


0.12% 


32.08% 


HOMOoc-5 


8.67% 


0 


8.08% 


0.07% 


0.11% 


43.12% 


HOMOoc-6 


9.96% 


0.01% 


15.20% 


0 


0 


0.06% 


HOMOoc-7 


25.27% 


0.12% 


29.89% 


0.03% 


0.01% 


2.89% 



orbitals localized on the U atoms has been deformed and no longer 
has any symmetry of U2@C6o^°. Furthermore, we can obtain the 
natural atomic orbitals (NAO) from the density matrix. These 
NAOs can be understood as deformed atomic orbitals in the molecu- 
lar structure. Molecular orbitals can be expanded by these NAOs, 
which are orthogonal to each other. So, we can calculate the sum of 
the square of the corresponding orbital coefficients as the contri- 
bution to these molecular orbitals. As the 5f and 6d atomic orbitals of 
U atom are not fully occupied, they usually have a significant con- 
tribution to bonding. The 2s and 2p orbitals are the valence orbitals 
of C atoms. These atomic orbitals are active in intermolecular inter- 
action. Table 4 shows the contributions in the form of a percentage. 

According to the shape of the MOs in Figure 4 and corresponding 
to the data in Table 4, all these MOs reflect the bonding between the 
two U atoms. The HOMOa-4 and HOMOa-5 reflect the contri- 
bution from the 2p electrons of the ad-atom. The 5f electrons of U 
atoms have a significant contribution to these frontier MOs but the 
contribution of 6d shell is very small. These results show that the 
eight singly-occupied electrons mainly come from the 5f of U atoms 
and the 2p of the ad-atom. 

We can directly find out the source of the contribution of the spin 
state from the spin density analysis. The ground state spin density 
analysis of U2@C6i (c, ^A') is shown in Figure 5. The ground state of 
the neutral U2 molecule is known to have a septet spin state. The six 
electrons of the same spin are contributed by the 5f orbital of the two 
U atoms. The ground state of C^i fullerene cages is triplet due to the 
existence of the defect, and the unpaired electron comes mainly from 
the 2p of the C adatom (spin = 1.82, supplementary information Part 
1.3). We can clearly see from Figure 5 that the nonet ground state of 
U2@C6i mainly gets the spin contribution from both the U atoms 
and the adatoms. The unpaired electrons on the U atoms and the C 
ad-atom are of the same spin (the color of them are all blue). 
Calculated results show that the spin density of Uist is 3.09, U2nd 
of 3.30 and the adatom of 1.51. And also, there is a small amount of 
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spin contribution arising from the other C atoms of the fullerene 
cage. So, the reason of the nonet spin state is that there is a ferromag- 
netic coupling between the spin of the unpaired electrons of the U 
atoms and the adatom. 

Discussion 

In summary, we have constituted an EMF of U2@C6i structure 
through a spin polarization adatom-type defect on the surface of 
the typical dimetalloendofullerene U2@C6o. DFT theoretical calcula- 
tions have shown that the defect could expand the local space within 
the Ceo fullerene. This expanded local space and the electronic 
impact of the ad-atom directionally 'attract' one of the two U atoms 
to the defect area. The bond length in U2 molecule elongates signifi- 
cantly in Cei fullerenes as compared to that of U2@C6o. From the 
vibrational spectra and molecular orbital analysis, bonding interac- 
tions are still observed between the two U atoms in C^i fullerene, 
signifying that the expansion of the local space in the Ceo cavity 
caused by the ad-atom defect does not dissociate the U2 molecule 
completely. Moreover, this EMF has a unique electronic state unlike 
the U2@C6o. The electronic ground state of U2@C6i has a nonet spin 
multiplicity with a ferromagnetic coupling between the spin of the 
unpaired electrons of the U atoms and the adatom. 

Because of its unique electronic structure, EMFs have very import- 
ant promising applications in functional nanomaterials, biomedicines, 
etc. Our study found that we can further modify the equilibrium 
structure and electronic structure of EMFs by the introduction of spin 
polarization defect. This work will contribute to a profound under- 
standing of the complex interplay between the inner metal atoms and 
outer surface of the fullerene cage in different defect-induced environ- 
ments or conditions. It can pave possible new ways to regulate the 
electronic structure of the EMFs to realize its practical applications. 

Methods 

The relativistic effects and electron correlation must be taken into the consideration 
in a systematic investigation of actinide- containing systems. A significant advantage 
of DFT calculations is that accounting for the electron correlation does not increase 
the CPU time considerably. Therefore, DFT has been widely used in the calculation of 
actinides^^"^^. Relativistic pseudopotential method is able to incorporate relativistic 
effects with a small increase in computational time^^'^^. The impact of the pseudo- 
potential core size of U on chemical calculations based on DFT has been reported^^'^*^. 
Previous studies show that pure DFT methods perform better than hybrid DFT or 
many-body perturbation methods in describing metal diatomics^^"*\ and therefore 
pure DFT, such as BP86 and PBE, was used in U2@C6o studies^^'^". Based on the above, 
the gradient-corrected BP86'^^ XC-functional was employed in this paper, and the 
hybrid functional PBEO*^'^^ have also been used for a quality check. Based on testing 
results and analysis of the structure and charge population characteristics, 3-2 IG 
basis set''^ was chosen for the carbon atoms^*^"^". Similarly, 3 -2 IG basis set for C atoms 
of fullerene has been applied in a previous lanthanide EMFs study^\ Scalar relativistic 
effects were applied to the uranium atoms. The (14sl3pl0d8f6g)/[6s6p5d4f3g] 
named ANO^^ valence basis set and ECP60MWB^^'^^ pseudopotential basis set with 
60 core-electrons were used. All the calculations have been performed in Gaussian 09 
program^^. 



Figure 5 | The spin density distributions of the U2@C6i (c, ^A'). 
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